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Abstract 

We obtain and analyze the indirect exchange interaction between two 
two-state systems, e.g., spins, in a formulation that also incorporates the 
quantum noise that they experience, due to a bosonic environment, for 
instance, phonons. We utilize a perturbative approach to obtain a quan- 
tum evolution equation for the two-spin dynamics. A non-perturbative 
approach is used to study the onset of the induced interaction, which is 
calculated exactly. We predict that for low enough temperatures the in- 
teraction is coherent over time scales sufficient to create entanglement, 
dominated by the zero-point quantum fluctuations of the environment. 
We identify the time scales for which the spins develop entanglement for 
various spatial separations. 
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Studies of open quantum systems have a long history ^E]. Recent exper- 
imental advances have allowed the observation of fundamental quantum me- 
chanical phenomena in nanostructured systems in condensed matter and other 
fields. Promise of applications for quantum information processing has stim- 
ulated significant interest in theoretical studies of quantum coherence and en- 
tanglement in situations when the quantum system is subject to environmental 
noise ^1 12 03 El 03 . I n the present work, we consider two two-state systems: 
qubits, e.g., spins 1/2, in a thermal bosonic environment (bath). We study 
the emergence of the indirect exchange interaction between two localized spins 
and identify the regimes where entanglement generated by this interaction can 
be observed. We report results for geometries relevant for recent experiments 

013111. 

Quantum computing schemes with qubits coupled by indirect exchange of 
excitons were proposed in OEHEJEL anc ^ these interactions are also of inter- 
est in studies of quantum phase transitions [13 . Traditionally, such RKKY-type 
interactions were calculated perturbatively at zero temperature, without con- 
sidering quantum noise, with conduction electrons ^3] or exitons ^B] acting as 
the "bath." Recent works ^3E| considered the effects of noise for a thermal- 
ized bath of, respectively, noninteracting and interacting electrons. A bath of 



thermalized bosonic modes can cause decoherence and, for more than one qubit, 
disentanglement. These effects have been studied extensively in the recent lit- 
erature It has also been anticipated p3] that such a thermalized bath of 
modes can induce entanglement under certain conditions. 

In this work, we investigate both physical effects of a thermalized bosonic 
bath in which two qubits are immersed. Specifically, the induced interaction, 
which is effectively a zero-temperature effect, and the quantum noise, originating 
from the same bath modes, are derived within a uniform treatment. We study 
the dependence of the induced interaction (coherent) vs. noise (decoherence) 
effects on the parameters of the bath modes, qubits, and their coupling, as well 
as on the geometry of the qubit system. 

We consider two localized spins separated by distance d and identically 
coupled with the modes of a thermalized bosonic bath, described by Hb = 

^kajcOk, where we set K = 1. The external magnetic field will be represented 
by the Hamiltonian Hs, corresponding to the energy gap A between the up and 
down states for each spin. A natural example of the described model are spins 
of two localized electrons interacting via lattice vibrations (phonons) by means 
of the spin-orbit interaction [201 • For each type of phonons, the interaction can 
be assumed ^ |3 E01 linear in the bosonic variables. Without loss of general- 
ity, the calculational techniques can be developed for a coupling that involves a 
specific spin component, which will be the same for both spins and denoted by 
Sj for the jth spin, localized at Tj. Thus we take Hsb — $Zj=i 2 SjXj, where 



Xj = X)k5ke ikri (ak + alkj- All our derivations, however, can be general- 
ized to include all the projections of the spin and we give some illustrative 
numbers for such a calculation below. Our emphasis here is on comparing the 
relative importance of the coherent vs. noise effects of a given bosonic bath in 
the two-qubit dynamics. We do not include possible other two-qubit interactions 
in such comparative calculation of dynamical quantities. 

The overall system is described by the Hamiltonian H = Hs + Hb + Hsb- 
Let us point out that such a model is quite general and it also finds applications, 
for instance, in quantum optics |22| where the Hamiltonian H would describe 
atoms (regarded as the two level systems) interacting with an electromagnetic 
field. Our detailed expressions here are obtained for the one-dimensional case, 
relevant for recent experiments that involve channel geometries [H1E1IH], when 
phonons propagate along d, i.e., k • d — * k\d\. 

The problem of describing the dynamics of the spin system and, thus, finding 
the reduced density matrix, ps(t) = Trsp(t), cannot be solved exactly in the 
general case. We first consider time scales longer that the thermalization times 
of the bath. Furthermore, we treat the interaction, Hsb, as a perturbation, and 
expand the equation for the density matrix, idtp(t) = [H, p(t)], up to the second 
order in Hsb- In this regime, it is appropriate [2j to model the thermalization of 
the bath within a standard Markovian approximation which involves factoring 
out the environmental mode density matrix in the second order term, replacing 
it by the thermal one, as well as using P the initial condition p(0) = ps (0) <g> 



. Here p"?™*' = rj fc z- 1 e~ UJka l ak / kBT , and Z fc is the partition function 
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for the oscillator k. The resulting equation [2] for the density matrix is 

id tPs (t) = [H S , ps(t)] + i I dt'V (f - t) Ps (t), (1) 

Jo 

where £ (f - t) Ps (t) = - Tr B [S 3 X 3 , [5< (f - t) X, (t> - i) , pf eTmal p s (t)] ] 
is the self-energy superoperator term. 

This expression involves the correlation functions Cji(t) — Trs[XjXi(t) 
XPs erm ], where j,j = 1,2 for the two spins, with the property C*,j(i) — 
Cji(—t). The present approximation assumes 21 that the bath has very short 
memory. One can argue that the correlation functions Cji(t) are nonnegligible 
only up to times 1/cj c , where uj c is the frequency cutoff for electron spins inter- 
acting with phonons. This cutoff comes either from the phonon density of states 
or from the localization of the electron wavefunctions. This time scale will be 
considered to be significantly smaller then the system evolution times defined 
by the inverse of the energy gap A. 

In the resulting evolution equation, 

id tPs (t) - [ffeff, Ps(t)] + i£cPs(t) + iC sPs (t), (2) 

we separate out the coherent dynamics in the first term. The superoperators 
C c and L s will be addressed shortly. In the two-spin case, one can establish by 
a lengthy calculation which is not reproduced here, that 

fleff = H s + 2 X c{d)S 1 S 2 + O ( Xs {d),r h (0)) , (3) 

which includes the effective coupling of the spin components. Expressions for the 
quantities \c,s & n d ?y c ,s follow from, respectively, the imaginary and real parts 
of the correlation functions and will be given explicitly below. The last term in 
P| introduces corrections of relative order A/lo c in the induced interaction and 
is not of interest here. 

The second term in J2J is 

£c P s{t) =-^Vc(9i j d){SiS J ,p 8 {t)}+2j^r h {e ij d)Sips{t)S i , (4) 

ij i^j 

where 6ij = 1 — 6^; it accounts for the dominant relaxation and decoherence 
processes. The third term in |J2j, C s ps{t), involves expressions proportional 
to r) s (Oijd) or Xs (@ijd) and may often be considered small compared to terms 
in J!}, as will be shown bellow by analyzing the magnitudes of rj s (Oijd) and 
Xs (Oijd). 

In and @ we introduced the quantities 

oo 

Xc (d) — J dt J du> V (u) \g (w)| 2 sin uot cos cos At (5) 

-oo o 
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and 



oo 

r\ c (d) = f dt ( dujV (to) \g (lu)\ 2 coth — - — — cos cot cos cos At, (6) 
J J 2k B T c s 

-oo 

where V (w) is the density of states of the bosonic modes; this factor should also 
include the Debye cutoff at large frequencies. For definiteness, we have assumed 
the linear dispersion, w = c s k, because details of the dispersion relation for larger 
frequencies usually have little effect on decoherence properties. Another reason 
to focus on the low-frequency modes is that an additional cutoff, ui c , resulting 
from the localization of the electron wave functions, typically much smaller than 
the Debye frequency, will be present due to the factors \g (w)| 2 . The expressions 
for rjs (d) and % s (d) can be obtained by replacing cos At — > sin At in r\ c (d) and 
Xc (d) , respectively. The integration in © then yields 

Vc (d) = |D (A) \g (A)| 2 coth ^ cos (7) 
Similarly, we find 

Xs (d) = ^( A )lf( A )| 2cos ^- (8) 

Z C s 

To derive explicit expressions for Xc (d), one needs to specify the w-dependence 
in V{u) \g(w)\. For purposes of modeling bosonic heat-bath effects, this prod- 
uct is usually approximated T by a power law with superimposed exponential 
cutoff, V (w) \g (w)| 2 = a„w"e-^, with n > 1. Finally, we get 



Xc (d) 



0"- 1 (-1)""^ 
5 d^- 1 l + t\ 2 



(9) 



This quantity gives the coefficient of the leading induced interaction; see ©• 
The individual Lamb shifts in @ are defined by r) s (0) and are not important 
for our discussion. 

For the commonly studied case of Ohmic dissipation PP, n = 1, the environ- 
ment induces the following spin-spin interaction, 

H m t — — ; 9 i9 / 9 ^1^2 ■ (10) 

f + bj^a 1 1 cj. 

This induced interaction is temperature independent and is mediated by the 
zero-point fluctuations of the bosonic field. It is long-range and decays as a 
power-law for large |d|. On the other hand, quantum noise terms, see II) . 
depend weakly on |d|, and increase with temperature. In Fig. 1, we plot the 
magnitude of the interaction Xc (d) as a function of the spin separation. It is 
compared, for varying temperature, to the magnitude of the decoherence terms, 
Xs (d), r/s (d) and r\ c (d), among which rj c (d) is the dominant for the Ohmic 
case, n = 1, and the temperature scale of Fig. 1, i.e. for 2kT/ui c from ~ 0.05. 
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For lower temperatures evaluation of rj s (d) may be important. For the super- 
Ohmic case, n > 1, one obtains similar behavior, except that the interaction 
decays as a higher negative power of |d|. 

The amplitudes plotted in Fig. 1 provide qualitative information on the 
dynamics. For defmiteness, let us consider the case of Hs — ^Aa^ + ±Aa^ 
and -ffi n t in (|1L)|I involving Si 2 = cri 1 ' 2 ^. If the induced interaction were the 
only effect of the bath, then the system's dynamics would be coherent, with 
oscillations determined by the energy gaps of Hs + -ffint- This Hamiltonian 
has the singlet state (|fj.) — ||T)) /V^, with the energy E2 = — 2% c (d), and the 
(split) "triplet" C(|tt )-*|U)), (IT 1) + |!T»/V2, and C(||l) + 5 |TT)), with 
the energies E = - y/ A 2 + 4 X 2 {d) , E x = 2 Xc (d), and E 3 = y/ A 2 + A\ 2 (d), 
respectively, where 5 = 2\ c (d) / (A — Eq) : and C is the normalization constant. 
The effect of the noise terms is to wash away the coherent behavior. The time 
scales of the coherent oscillations and of the noise-induced relaxation processes, 
defined by the corresponding frequencies in Fig. 1, become comparable at the 
intersection curve of the two surfaces in Fig. 1. For larger distances and/or 
temperatures the noise terms dominate. 

The energy gap between E\ and E 2 is defined by Xc (d), whereas the effective 
width of each level due to decoherence will be determined by the magnitudes of 
Xs (d), rj s (d) and r\ c (d). Let us estimate the magnitude of this energy splitting. 
As an example, we consider two phosphorus donor impurities in Ge, in external 
magnetic field in the z direction. Due to the symmetry of the spin-phonon 
coupling (via spin-orbit interaction), only two components, namely, a x and a v , 
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are important [53]. One can demonstrate that in this case the terms in |QJ 
proportional to products of a x and a y vanish, while the contributions from the 
<j x and <7 y terms are identical. Therefore, we can use the single-component 

results, with Sj ~ Ox . For a donor impurity electron spin, the cutoff oj c ~ 
c s /cib comes from the donor electron wave function |23j . which is localized on 
the scale of clb ~ 4nm. The characteristic value of the phonon group velocity 
in Ge is about 3 x 10 3 m/s. The strength of the spin-orbit coupling, aiw c , is 
of the order of 10's -1 . Utilizing these data in (|7I10|I . with n = 1 which, 
strictly speaking should be only valid for one-dimensional channel for phonon 
propagation, we obtain the estimate for the E\ <-> E 2 energy splitting of about 
10 MHz, whereas the noise level for mK temperatures varies from 0.1 to 1 MHz, 
depending on the magnitude of the Zeeman splitting. This coherence/noise 
"measure" can be further improved for different cases of the phonon spectrum, 
the shape of the wave function, etc. |21) . 

To study the onset of the exchange interaction, we consider the case when 
the Zeeman splitting A is negligible. Then one can actually derive an exact 
solution for H$ + Hg + Hsb, with the bath modes traced over without the 
Markovian assumption, and demonstrate the emergence of the effective inter- 
action Hamiltonian as in I|3I9I10|) . For A = 0, a lengthy calculation utilizating 
bosonic operator techniques yields 

Ps(t) = Y,P\Ps(0)Pye c ^' (t \ (11) 

A, A' 

where |A,-) is an eigenstate of Sj labeled by its eigenvalue A^, and we introduced 
the projection operator P x = | A x A 2 ) (Ai A 2 1 ■ The exponent consists of the real 
part, which represents decoherence, 

ReC xy (t)= - Y, G ^ T ) [(Ai - Ax) 2 + (A' 2 - A 2 ) 2 
fe 

+ 2cos(w fc |d|/c.)(Ai-Ai)(A£-A 2 )], (12) 
and imaginary part, which describes the coherent evolution, 

Im£ AV (t) =J20 k (t)co S (uj k \d\/c s ) (AiA 2 - X[X' 2 ) . (13) 

k 

Since for A = 0, Hs commutes with Hsb, quantum noise in this case only 
affects the off-diagonal matrix elements. Eventually, this destroys quantum 
correlations between the spins. In l(T2l and the standard functions PHJ|3] 
were introduced, 

G k (t; T) = 2^ sin 2 ^ coth (uj k /2k B T) , (14) 

Cfc(i) = 2-^- (uj k t - s'muj k t) . (15) 
u k 
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Figure 2: The onset of the indirect spin-spin interaction is measured by the 
decay of the correction term Fit). The values of F(t) are in units of a\uj Cl and 
are color-coded according to the top bar. 
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Figure 3: The concurrence calculated for Ohmic dissipation with a.\ = 
ksT/tOc = 1/20. The inset gives the time dependence for different tempera- 
tures: 80kBT/uj c = 1,2,3,4,5,6,7,8, from the top curve to the bottom one, 
respectively: the top curve corresponds to the lowest temperature. 



Our focus here is on the imaginary part (|13J) . One can demonstrate that if 
the noise terms (the real part) were completely absent, the resulting evolution 
would be coherent with the evolution operator given by e-^t+n*)]*, where 



F(t) = 25i5 2 



, V (no) \g(u>)\ smut 
duj — n — cos(w|d|/c s 



and 



2a n T (n)h£ 
(l+w»dVc2 



i/2 



tot 



n arctan 



S1S2 



(16) 



(17) 



Expression l|16|) represents the initial, time-dependent correction present only 
during the onset of the induced interaction. Specifically, -F'(O) = —H ln t, but 
F(t) decays for large times. In Fig. 2, we plot F{t) for the Ohmic case, n = 1, 
as a function of time and spin separation. The right side of the plot in Fig. 2 
corresponds to the constant coupling regime with the interaction Hamiltonian 
l|17|) . Note that the Hamiltonian l|17|) with n > 1 is identical to that obtained 
within the perturbative Markovian approach, cf. <|3I9I10|) . 

Let us investigate the role of the decoherence resulting from the real part of 
C\\> in (|12f) . Note that within the exact solution the bath is thermalized only 
initially, at t = 0. However, it is expected that the effects of the quantum noise 
are represented qualitatively similarly to the Markovian approximation, sec 

which implies rc-thermalization of the bath after each infinitesimal time step. 
To analyze the dynamics of quantum correlations between the spins, we have 
used the concurrence |24| — a measure of entanglement which is widely used 



in quantum information theory. For a mixed state of two qubits it is defined 

123] via the eigenvalues Xj of \J \f~Ps (J y ~y P*s a l a y\fPS as max{ 0, 2 max Aj — Ai — 

A2 — A 3 — A 4 }. For A = 0, and as long as the effects of the quantum noise are 
small, the induced interaction will split the system energies into two degenerate 
pairs, Eq — E2 and E\ — E%. As a result, the dynamics will involve coherent 
oscillations with the frequency defined by the gap E\ — E — 4% c (d). In Fig. 3, 
we plot the concurrence for the density matrix given by as a function of 
time and spin-spin separation, for the (initially unentangled) state ||T), and 
n = 1. The inset in Fig. 3 shows the time dependence of the concurrence for 
different temperatures. Figure 3 demonstrates that the system can develop and 
maintain entanglement over several coherent-dynamics oscillations before the 
noise-induced effects take over and the concurrence decays to zero. 

In summary, we studied the induced indirect exchange interaction due to 
a bosonic bath of environmental modes which also introduce the noise. We 
demonstrated that it can create observable two-spin entanglement. For an ap- 
propriate choice of the system parameters, specifically, the spin-spin separation, 
this entanglement can be maintained and the system can evolve approximately 
coherently for many cycles of its internal dynamics. However, for large times 
the quantum noise effects will eventually dominate and the entanglement will 
be erased. 

The authors acknowledge useful discussions with and helpful suggestions by 
J. Eberly, L. Fedichkin and D. Mozyrsky. This research was supported by the 
NSF under grant DMR-0121146. 
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